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Abstract. We calculate the trispectrum of curvature perturbations for a model of 
inflation endowed with Galilean symmetry at the level of the fluctuations around an 
FRW background. Such a model has been shown to posses desirable properties such as 
unitarity (up to a certain scale) and non-renormalization of the leading operators, all 
of which point towards the reasonable assumption that a full theory whose fluctuations 
reproduce the one here might exist as well as be stable and predictive. 
The cubic curvature fluctuations of this model produce quite distinct signatures at 
the level of the bispectrum. Our analysis shows how this holds true at higher order 
in perturbations. We provide a detailed study of the trispectrum shape-functions 
in different configurations and a comparison with existent literature. Most notably, 
predictions markedly differ from their P(X, <j>) counterpart in the so called equilateral 
trispectrum configuration. 

The zoo of inflationary models characterized by somewhat distinctive predictions for 
higher order correlators is already quite populated; what makes this model more 
compelling resides in the above mentioned stability properties. 
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1. Introduction 

The inflationary paradigm [1] stands at the core of modern cosmology. It set out to 
solve several cosmological puzzles which pre-dated it and delivered further detailed pre- 
dictions verified to be in close agreement with observations [21 EJ [U [5] . The simplest 
models of single-field slow-roll inflation are, as of today, in agreement with data. On the 
other hand, it is natural to ask whether it is possible to find or embed an inflationary 
mechanism within UV complete theories such as string theory. One such a realization 
(DBI inflation) was put forward in [HI IB] , and has since been studied in great details 
(for other realizations see [9], jTOJ HH [T21 [13] and [THEE! [16] f° r some reviews). However, 
the regime in which the DBI model can be trusted seems perilously close to what data 
might soon exclude. 

Stopping a bit short of this last, quite ambitious, requirement, one could also aim at 
a model of inflation which is unitary, stable under quantum corrections and predictive. 
In [T7] the authors propose an inflationary set up in this direction, where the require- 
ment of the so called Galileon symmetry (0 — > <p + c + b^ , see [19] ) grants all of the 
above properties: 

- in flat space the n (where the index n counts space-time dimensions) Galilean terms 
are guaranteed to have second order equations of motion [12] 

- the non-renormalization theorem for Galileon models [19J warrants that (at least in 
flate space) the coefficients which one writes down at first will not be largely modified 
by renormalization. 
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Crucially, these desirable properties are shown to approximately (up to A/Mpi 
corrections, where A is the scale of the underlying theory) hold for the model in [17] 
also in the more realistic scenario where one deals with curved spacetime and couples 
to gravity. There exists an interesting energy regime where the Galileon interactions 
are important and the model is quite different from canonical slow-roll inflation. In 
this same regime both the metric fluctuations and terms with non-minimal coupling to 
gravity can be neglected and calculating higher order curvature correlators is greatly 
simplified. 

The bispectrum of this model has been explored in [17] : the non-linearity parameter 
/jvl generated can be detectable but the shape- functions peak in the equilateral 
configuration and thus one is not able to remove the large degeneracy with most other 
inflationary models at the level of the bispectrum. This, of course, leaves the appeal 
of the model just described untouched. One should then go look for another quantity 
at (possibly) observational reach, the trispectrum being the natural choice. This is 
precisely what we do in a paper [20J which is complementary to the work presented 
here. 

On the other hand, if one is prepared to give away a little bit in terms of the cher- 
ished properties described above (this statement will be made more quantitative below), 
then there exist a model [21] which already at the bispectrum level generates interesting 
and quite distinct shape signatures. Combining the analysis of [21] with the study of 
the four point function, the trispectrum, one can really hope to fully characterize the 
predictions of this generalized Galileon model and disentangle it from a host of infla- 
tionary models with different, somewhat more canonical, signatures. 

Briefly, in going from the approach of [IT] (see also [18]) to the one of [21] one: 

- has a smaller scale A up to which the effective theory can be trusted 

- with some assumptions (reasonable in light of [IT] and other works), preserves the 
non-renormalization properties of the Galileon terms 

- has a richer and more distinct set of predictions for non-Gaussianities 

What has been proven concerning the background (and therefore the full theory) in [TTJ 
is more or less implicitly assumed at the onset of [21] because the latter is an effective 
field theory of fluctuations around a background rather than a effective field theory as 
a whole. The focus is on fluctuations and their properties. 

Following [21] then, we study the four-point function of curvature perturbations in 
detail, we calculate the corresponding amplitude and plot the shape-functions of each 
interaction term in different configurations so as to ease the comparison with the existent 
literature on the subject and, at least at the level of shape-functions, remove as much 
degeneracy as possible. 

Our analysis nicely fits within the realm of the effective field theory approach to 
inflation |23 [2HH2HH21lI3ni[3IlE2l[331I31lES]- See [21] for a different approach. 

On the other hand, the model we study is not quite as generic as it could have been 
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in the spirit of [26J: we are now dealing with a theory which is unitary (up to some 
scale) and has been further endowed with a symmetry which selects interaction opera- 
tors and protects them from quantum corrections. This subset of theories is more stable 
and more predictive in that the symmetry shields us from large renormalization of the 
coefficients of interaction terms and keeps the number of leading interaction operators 
to a finite, small number. 

The paper is organized as follows. In Sec. 2 we introduce the action for the model 
and expand it up to fourth order in fluctuations, we also perform a field redefinition 
which returns the action in a form that makes it suitable for employment within the 
Hamiltonian formulation of the Schwinger-Keldysh formalism. In Sec. 3, we compute 
the trispectrum of curvature fluctuations, study its amplitude and shape functions so 
as to draw a comparison with some of the results already present in the literature for 
other inflationary models, especially DBI inflation (for studies on primordial trispectrum 
in DBI, and more general P(X,4>) models, and other inflation models, see, e.g., 
[391 SOI S2l H3l SSI S5l SZl SHI SH])- In Sec. 4 we conclude, and offer some 
further comments. In Appendix A we present the details of our field redefinition 
while in Appendix B we review the main steps that characterize the IN-IN formalism 
computation of the trispectrum and we list the final analytical results for the shape 
functions. In the following greek indices run from to 3, latin indices are spatial only, 
running from 1 to 3. 



2. The Action to fourth order in perturbations 

2.1. Building Blocks of the Action 

We are interested in a Lagrangian for scalar perturbations ir around an FRW background 
whose fluctuations satisfy the Galilean symmetry. The model can be written down by 
listing all possible Galilean invariant terms which non-linearly realize Lorentz symmetry. 

Alternatively, one can write down Lorentz invariant contributions in terms of 
ip = t + 7r(x, t) where ip linearly realizes the symmetry [21]. We shall proceed according 
to the latter prescription. What we are after is the number of independent operators 
generated by products of traces of the form: 

[^fci x ^jfa x ^ 

where each [...] is to be understood as a trace of whatever product is inside the brackets, 
^ = V^V^V) an d indices are raised with the metric g^ u . It can be shown that, once all 
the terms of the form [\l/ n ] with n < 4 are known, all the rest in Eq. ([T| can be expressed 
as linear combinations of products of these building blocks [2T] ■ 

With fluctuations in mind, one removes the background value c n from each trace below: 
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# 2 



- c% = V„VV - ci = V^W 
-C2 = V^ttV^tt - 2HVjW 
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a 3 
a i 



\[/ 3 



- c 3 = V^ttV^ttV^tt - 3#V Mi 7rV iA V + 3# 2 V;VV 



c 4 = ... - ^V^TrV^TrV^Tr + 4H 2 V ^vrV^vr + 2H 2 V lJ TxV' ji ^ - 4# 3 V»VV 

(2) 

where the expressions of the type V MlAt2 7r ; V^7r are meant as shortcuts for 
V^V^ff ; V^V^ 12 ^ (greek indices run from to 3, latin indices are spatial only) . 
It is easy to check that Eq. (J5]) represents a basis which spans the entire space (up to 
order 7r 3 ) of the [\l/ n ] — c n operators for any n > 4 . 

In what follows we do not consider contributions coming from a single trace e.g. [\I/ n ] ; 
the reason is that these will account |21j for the so-called minimal covariant Galileon 
terms [221 1231 [21] and those, their expression being known, can be put in at a later 
stage. More explicitly, the single trace operators can be written as sums of multiple 
trace operators and minimal Galileon terms. The former are already all accounted for 
in Eq. ^ while the latter can be shown to be subdominant whenever generic (V 2 ip) n 
operators are switched on, as is the case here [2T] . 

With Eq. g handy, one sets about to write down all the leading independent 
interactions at third and fourth order in perturbations. Particular care must be 
exerted so as to have that, at each perturbative order, the coefficient multiplying a 
given interaction term does not appear already at lower perturbative order. Such an 
occurrence is not problematic by itself but, as is known [2Tj . the third and fourth order 
contribution of a piece like M 4 (n 2 + 7r 3 + 7r 4 ) would be suppressed by a scale A' larger 
than the one A < A' that suppresses the quadratic fluctuations. The resulting three 
and four-point functions would therefore be negligible with respect to the contribution 
coming from interactions suppressed by A. 
Indeed, what is natural to assume at this stage is that: 

Every independent interaction term at each given perturbative order is suppressed by 
a common scale A . This is possible if each coefficient which corresponds to each 
independent interaction appears at given perturbative order and nowhere else [|] 

This requirement is relatively easy to implement. Using the building blocks in 
Eq. g, what one wants is independent combinations of product of multiple trace 
operators whose first contribution starts at a given perturbative order 

= E Cm(x!,x 2 , x 3 ) x a^afafa^ = M x 0(ir 3 ) + ... (3) 

X 

04 = E Cm(xi,x 2 , x 3 , X4) x al 1 a x 2 2 al 3 al 4 - = M x 0(tt 4 ) + (4) 

X 

where the c's are constant (to first approximation in generalized slow-roll) coefficients. 
As an example, consider: 

a\ ~ M 4 (ViVV) 3 ; a\ ~ M 5 (V^7r) 4 ; (a A + \Ha^) ~ Mi (v^-vrVV) 2 . (5) 

The coefficients M&, M n are all different and appear only once in our total Lagrangian. 
Applying this algorithm to the fullest, one obtains the interactions Lagrangians below. 

| That is, appears as multiplying a leading operator. 
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2.2. Interaction Lagrangian 

We write below the expression for the 1S3, £4 action in tc: 

«S 3 = / dtd'x^-giM^ln) 3 + M 2 (V 2 7r) 2 (V 2 7r) + M 3 (V 2 7r)(V 2 7r) 2 + M 4 (V 2 tt) 3 

+M 5 [(V M ,7r) 2 - 2(V p ,7r) 2 ](V 2 7r) + Af 6 [(V^7r) 2 - 2(V p ,7r) 2 ](V 2 7r) 
+M 7 (V u -vr) 2 (V 2 7r) + M 8 (Vy7r) 2 (Vfr)} (6) 

S 4 = J dtd 3 x^g~{M 1 [(V ij n) 2 } 2 + M 2 (Vy7r) 2 (V 2 t 7r) 2 + M 3 (VJtt) 4 + M 4 (V 2 vr) 4 

+M 5 (V 2 vr) 3 (Vjvr) + M 6 (V 2 7r) 2 (V 2 7r) 2 + M 7 (V 2 7r)(V 2 7r) 3 
+M 8 (V 2 7r) 2 (V, fc vr) 2 + M 9 (V 2 vr)(V, fc vr) 2 (Vj7r) 

,2 / , ,„\2 
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((V M ,7r) 2 ) + 4 ((V^vr) 2 ) -4(V, p 7r) 2 (V /ll 7r; 



+M 11 [(V^ - 2(V /9J tt) 2 ](V» 2 + M 12 [(V^tt) 2 - 2(V pi 7r) 2 ](V^ 
+M 13 [(V At ,vr) 2 - 2(V pi 7r) 2 ](V 2 7r)(V 2 7r) + M 14 [(V^vr) 2 - 2(V p4 vr) 2 ](V jfc 7r) 2 } 



(7) 



The second order Lagrangian is 

S 2 = f dtd 3 x^gM 2 P H (^vrd'V) . (8) 



If one were to write down the ir equations of motion for, say, S3 above, those would 
surely include terms like third time derivatives of the scalar d 3 ir/dt 3 ] that is because 
the action above comprises higher time derivative terms. If taken literally, this fact 
leads in turn to the excitation of unstable modes in the 7r dynamics and renders the 
Hamiltonian of the system unstable [SU]. In perturbation theory, it is however possible 
to treat this instability by essentially connecting the "higher derivative expansion" with 
the perturbative expansion in the field (52]. In the case at hand the latter procedure 
corresponds to performing an order-by-order (in perturbations) field redefinition. If one 
stops at third order in fluctuations, it can further be shown that a field redefinition is 
equivalent (up to C(vr 4 ) corrections) to the repeated use of the iS 2 -derived equation of 
motion to get rid of the time derivatives in excess. This is where we part ways with 
the analysis in |21j: there the bispectrum was the target and the use of the e.o.m.'s 
was entirely justified to this aim, here we also want to calculate the trispectrum so we 
have to undergo (as also recognized in [21]) a full fledged field redefinition: the C(vr 4 ) 
corrections do matter for our purposes. 



Let us redefine the field ir as follows 

7T ->> 7T( ) + 7T(i) + 7T( 2 ) + (9) 

where 7T(i) ~ O (flfo))) 7r (2) ~ O (flfo))- We are computing the action up to fourth order 
in the field fluctuations, therefore we do not need to go any further than 71Y2). If we plug 
the expansion in Eq. ([9l in the action up to fourth order in 71", we have 
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S 2 [n] S?[nf 0) ] + S?[n (0) n {1) ] + S?[7r 2 (1) ) + S? ] [n m n (2) ] + (10) 
S 3 [n]^S^f 0) ]+Sf\rr 2 0) 7r {1) ] + ..., (11) 

S 4 [7r]^4 4 H)] + ---, (12) 
where, on the right-hand side, the lower index on S indicates the perturbative order 
whereas the upper index reminds us from which order in S[tt] that specific contribution 
arises. The indicates contributions that are beyond fourth order and which we do 
not need for our computation of the trispectrum. The third and fourth order action in 
the newly redefined field 7T(o) are then 

s 2 [7r ( o)] = 4 21 [4)]> ( 13 ) 

«S 3 [7r (0 )] = «S| Jl [7rf 0) ]+5f[7r ( o)7r (1) ], (14) 

<s 4 [7r (0 )] = <si 41 [4 } ] + 4 2 W)] + s?W (a) ] + (is) 

We will now determine 7T(i) and 7T(2) in such a way that the final equations of motion will 
be at most second in time derivatives and, in the process, derive the final expressions 
for the third and fourth order Lagrangian. 

2.2.1. Cubic Lagrangian 

Upon expanding Eq. ([6]), one can verify that the terms in iS 3 that are expected to give 
rise to unwanted higher order time derivatives in the equation of motion for 7r are 

4 3 V( 0) ] 3 / d 4 xa 3 {nf 0) (-M x - M 5 )+n 2 0) (^^-3Hn {0) ) (ZM X + M 2 + M 5 + M 6 )}.(16) 

It is straightforward to check that the following field redefinition grants a cancellation 
of the above terms thanks to some of the contributions from 5 3 2 ' [^{0)^(1)} 

'Mi + M 5 \.. 2 (3H (2Mi + M 2 + M 6 ) \ (<IM X + M 2 + M 6 \ .. d 2 n 



(o) 



71 w = [ ^pn, ) ^(0) + 1 1 nmo) - [ — — ) n (°y „ 

(17) 

The final expression for the action to third order is obtained by summing the leftover 



terms from the two contributions in Eq.(14) (while the fourth order Lagrangian is 



generally affected by 7Tm, no contributions to the cubic Lagrangian can possibly arise 
from 7T( 2 )) 

SaM = / rf 4 xa 3 {Q^ (0) - 3Hir {0 ^ + Q 2 - 3Hn (0 ^j 



where 

Qt = -Mi -M 2 -M 3 + M 6 , (19) 

()■, ■ M 2 • .1/,, • M 4 , (20) 

Q 3 = M 5 - M 7 , (21) 

Q 4 = -M s - M 6 + M 7 + M 8 . (22) 
We can now move on to the quartic Lagrangian. 

2.2.2. Quartic Lagrangian 



The fourth order action is computed from Eq. (15) with 7T(i) as given in Eq. (17). The 
expression for km) is to be determined in such a way that some of the contributions from 
<S^ [ 7r (o) 7r (2)] will cancel the higher order time derivative terms from <S| [7tf J, S] 2 [^m] 

[3]r 



and <Si J [7rf Q) 7r ( i)]. 

The process of computing tt^) is straightforward but rather long, therefore we report 
the details in Appendix A. Before providing the final result for the quartic Lagrangian, 
it is useful to make some considerations about the mass scales characterizing the 
action. Let us consider our initial Lagrangian (j^-Q, with coefficients M and M, 
respectively for third and fourth order contributions. The field redefinition (in the 
form of 7T(2)) introduces extra interactions at the quartic level that are proportional to 
M = (M 2 H 2 )/(MpH). Therefore, our Lagrangian after field redefinition has the form 

C ~ M 2 H (dit) 2 + M (<9 2 7r) 3 + M {d 2 ixf + M (<9 2 vr) 4 . (23) 
Introducing the canonically normalized field ti c = MpH 1 ' 2 ^, we have 

„ , n ^2 M / 9 X3 M / n9 \4 M 2 H 2 / n , a . . 

C ~ (d-n c f + —. — - (d 2 7r c ) + (d 2 7r c ) + s (d 2 7r c ) ■ 24 

V ; (M 2 Hf/ 2 ^ > M P H 2 ^ > (m 2 E) K ' 

We restate the assumption (as in [21] ) that all operators multiplied by a coefficient that 



appears at a given order for the first time, are suppressed by a common scale A. Eq. (24) 
can then be rewritten as 

c „ (ftrJ » + ^ (a\^ s + 1 (a^f + Aj (av) 4 , (25) 

where we defined the scales 

k ^M i H^ = MFHV 1 vJj&V. (26) 

M M 5 / 8 ' M 2 H 2 y 1 

At this stage one can quickly relate the expression for M to the one for M 



A 
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M = M x ( — J M > M, (27) 

8 



where A ^> H if we are to make sense of perturbation theory. This is equivalent to 
stating that the scales of suppression are related by A > A. In other words, at each 
perturbative order, all operators multiplied by coefficients which appear for the first 
time at that order, are leading compared to operators that are multiplied by coefficients 
that appear for the first time at a lower order. Based on these considerations, we can 
safely neglect all operators proportional to M in the final fourth order Lagrangian which 
then reads: 



J a 8 a 8 a° a D 



r2 • 2 



V Q&n 7r - h Qjn tx — r- 



a" 1 a 11 a z 

+Q 8 ( -^P^ + Q 9 HW}, (28) 
a 8 

where Qi,---,09 are linear combinations of the initial coefficients Mj (we report their 
exact definitions in Appendix A). The final expression for the quartic action is now free 
of higher order time derivative terms and suitable for utilizing the Schwinger-Keldysh 
formalism j36j E3 EE] in the standard Hamiltonian formulation (as in [ID]) for the 
computation of the trispectrum of curvature fluctuations. Eq. (28) is one of our main 
results. 

We want to stress here that the field redefinition has not rendered the theory ghost- 
free, our procedure simply organizes the perturbative expansion in such a way as to have 
perturbations (up to a given order) to be consistent up to the scale A. This is, after 
all, just an effective theory; if one wants to take the theory beyond this point, that is, 
above A, one should first provide arguments as to why the ghost is not there also in 
that regime. 



2.2.3. Interactions Hamiltonian 



At third order in the field fluctuations, the interaction Hamiltonian is given by T-Lmt = 
—Hint- This relation does not generally hold true at fourth order. We quickly review 
the procedure for deriving "H^ (see e.g. [ID], whose notation we borrow in this section). 
Consider a Lagrangian of the type: 

£ = M 2 + j 2 + g ir 3 + gin 2 + g 2 n + Js + h n A + hn 3 + h 2 n 2 + h 3 n + j 4 , (29) 
where fo, gi, hj, are function of tt and its spatial derivatives (the subscript indicates 
whether a given function is order zero (/ , •••) linear (gx, ...) and so on in n). Introducing 
the conjugate momentum 

Ps w (30) 

we can derive the perturbative expansion, up to third order, of % in terms of powers of 
P. The total Hamiltonian reads: 

H(tt,P) = P-k - C, (31) 
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which is split into quadratic and interactions parts 

H = Ho + 'Hint- (32) 

In the Schwinger-Keldysh formalism, Hint is a function of tc and P as in the interaction 
picture ( tt 1 and P 1 ). Replacing 

71 =ap7> ( 33 ) 
in the interaction Hamiltonian, one can derive the following relations 

Hint = ~£-\nt-> (34) 

= "42 + O ( 9 -f). (35) 



/o 



In our specific case, the corrections to —C^ t in (35) are proportional to M — 
(M 2 H 2 ) / (MpH) . From our previous discussion on mass hierarchies, which lead to 
the relation M 3> M, we conclude that these corrections are subleading compared the 
terms in C^ t (regulated by M). In what follows we will therefore set Hi nt = —C int for 
the quartic order interaction Hamiltonian. 



3. Trispectrum of curvature fluctuations 

We now have all the ingredients necessary in order to compute the trispectrum of 
curvature fluctuations 

(kkkk) = (2vr) 3 ^ 3 )(fc 1 + k 2 + k 3 + h)T c (h, k 2 , k 3 , k 4 ) (36) 
where we Fourier expanded the curvature fluctuation (: 

(2^F e ^ (t) - (37) 

Whenever metric fluctuations are neglected, as is the case here, the curvature ( is linearly 
related to our field fluctuations n 



C = -Hn. (38) 

In particular, at leading order approximation in slow-roll, one can then essentially read 
off the ( interactions from the n Lagrangian. Upon quantization, the curvature is 
expressed in terms of creation and annihilation operators 

4 = a fc < fe + a^Cfc- (39) 

Using conformal time variable r, it is easy to derive from the quadratic fluctuations that 
the eigenfunctions are 

Ur) = - H ^^^ k \ (40) 

where Bunch-Davies vacuum has been assumed. The power spectrum of ( reads 

H 2 

AMlk 3 e 



(44> = (2«)*5®(k + k')P ( (k) => P C = (41) 
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(3) 

Figure 1. Diagrammatic representation of trispectrum contributions from H int 
(scalar-exchange, left figure) and H\„ t (contact-interaction, right one). 

The diagrammatic computation in the Schwinger-Keldysh formalism is outlined in 
Appendix B. At tree- level, is given by the sum of two distinct diagrams: a so-called 
scalar-exchange diagram, from ?4nL and a contact-interaction diagram from ?4nt ( see 
Fig. 1). It turns out that the former contributions are subdominant compared to the 
latter. This can be easily shown by comparing the trispectrum amplitude from the two 
diagrams. 

The amplitude of the trispectrum is conventionally described by a parameter, t^l, 
defined as the result of the normalization of the trispectrum to the third power of the 
power spectrum. r NL has a specific momentum dependence (related to the particular 
form of the interaction Hamiltonian) that we can only find by computing the interaction 
diagrams. 

However, knowing the expressions of the eigenfunctions (k and of the interaction 
Hamiltonian, dimensional analysis quickly leads us to the parametric dependence of 
tnl, which is necessary in order to estimate how big a contribution to the trispectrum 
a given diagram will provide. If we identify the contributions to the amplitude from the 
contact-interaction and from the scalar-exchange diagrams respectively as r NL and t^ l , 
we have 

H 2 M se H 2 M 2 

r NL ~ r NL ~ (42) 

Let us compare t nl and rff L : 

MpH 



which we have shown to be true (Eq. (27)). The contribution to the total trispectrum 
coming from terms represented by the contact-interaction diagram is then much larger 
than the contribution from the scalar-exchange side. We will therefore focus our 
attention on the former type of terms in the rest of the paper. 

3.1. Shapes and Amplitudes from contact-interaction diagrams 



Let us rewrite the Lagrangian in Eq. (28) in a schematic form 



11 



Sa 



d 4 xa 3 



(44) 



where we indicated operators by Oj (e.g. Q% = (dfn) 4 /(a 4 ) and so on). We calculated 
the contributions from all of the Oj interactions . Their analytic expressions are provided 
in Appendix B in the form of the shape functions TV * 
The the total trispectrum then reads: 

9 



see Eqs. (115) through (122)) 



■e, 
C • 



(45) 



The Tf *'s depend on the external momenta fci,...,fc 4 



which, because of momentum 



conservation, are oriented so as to form a closed tetrahedron. The total number of 
independent variables is six. We select them to be k%, k 2 , k 3 , fc 4 , k\2 and fcu, where 
ki = \ki\ and kij = \ki+kj\. One of the ki may be factored out in an overall normalization 
factor (we are in a scale invariant case), leaving us with five variables. This is of course 
too large a number for a plot of the momentum dependence of the 

One possibility is to select several different configurations in which the number of 
variables is narrowed down to two. We consider the configurations defined for the first 
time in |47j and applied to P(X, <fi) models therein (and later often times employed by 
other authors, see e.g. [29J [53]), so as to provide a clear comparison between results 
for the trispectrum in our model and the corresponding outcome in other inflationary 
models. 



A first glance at our results, from the qualitative appearance of the plots of the T®* 
for each configuration, one realizes the shape functions can be classified in two main 
categories. The analytical significance of these two categories is related to whether the 
respective vertices include or not the operator (didjir) 2 . It is also important to observe 



that the coefficients Q\ through Q 9 in Eq. (28) are not all independent from one another 



In fact, as one can easily check by looking at their definitions in terms of the original 
coefficients Mj, only four out of the nine Qi are independent. One therefore expects 
no more than four independent contributions to t^l- We select a basis of independent 
coefficients and collect the operators accordingly. 

The coefficients {Qi, Q2, Q7} are linearly independent from one another, therefore 
they represent a natural candidate basis for the ensemble Qi, ...,Qg. In this basis, the 
Lagrangian becomes 



where 



Sa — I d 4 xa 3 



Q 1 e 1 + Q 2 e 2 + Q 4 e 4 + Q 7 e 7 



0i = 0i- 12 3 + 54 5 - 162 



9- 



2 EE 2 - 2 3 + 9 5 - - 8 - 27 9 , 



(46) 

(47) 
(48) 



12 



e 4 = e 4 - e. 



1 9 

4®* + a 6 * 



9 7 ee9 7 



9 

7 9- 

4 



(49) 
(50) 



From the analytic form of the diagrams with vertices G7 and 0g (Eq. (121)), we can 
immediately notice that O7 gives a null contribution to the trispectrum. This means 
that we are left with three independent coefficients instead of four, which will generate 
three distinct contributions to t^l, along with the corresponding shape functions. These 
shape functions are plotted in Figs. 2-3. 
The expression of the total trispectrum is now: 

T c (k u k 2 , k 3 , h) = Tf 1 + Tf 2 + 2f 4 . (51) 

The form of each of the three contributions is 
T§1 = ( H*QA Tj 



(52) 



where the J-^s are functions of momentum which can be read off from Eqs. ( 115 ) through 
(122), using Eqs. (47)- (49). We consider the four momentum configurations introduced 



in [47] . It turns out that, in all configurations, T 2 and J-4 are qualitatively very similar 
(although their analytic expressions do not coincide), whereas T\ differs from them. 
As anticipated before Eq. (46), these similarities/differences are related to the operator 
(didjir) 2 : the vertex producing the T\ contribution, unlike the ones producing T 2 and 
J-4, does not contain (didjTt) 2 . We will comment on the plots of T\ and P 2 , in Figs. 2 
and 3 respectively (all the features of P 2 are shared by J-4). 
We provide below a configuration-by-configuration analysis of the results. 



Equilateral configuration. The four external momenta have equal length {k\ = k 2 = 
^3 — ^4)- The leftover variables are k\ 2 and ku, which we normalize to k\. Triangular 
inequalities restrict their domains. If we define x = k\ 2 jk\ and y = ku/ki, then 
x,y G [0,2] and y < y/A — x 2 . As we can see from Fig. 1, T\ is constant in this 
configuration, a functional behaviour that is also typical of contact-interaction diagrams 
in P{X,<j>) models g7]. 

On the other hand, the shape of T 2 (Fig. 2) is not a plateau and it very much 
resembles a shape that arises from contact interaction diagrams with the operator (V7r) 4 
in Ghost Inflation [53J; a similar shape was found within the Effective Field Theory 
(EFT) approach to single-field inflation in [29J, both from scalar-exchange (arising 
specifically from -h (didjir) 2 interactions) and from contact-interaction contributions (e.g. 
from (<9;7r) 4 and (d 2 n)(didj7rdjd k Ti d k diix)). 

We stress here that the aim of the paper is manifold. One is certainly after distinct 
features in higher order correlators, and indeed, we do find here signatures which are 
not present in quite general models such as P(X, 0) models. On the other hand, we 
complement these findings with the fact that they occur within a theory which is, from 
the effective quantum field theory point of view, stable and predictive. In some respects 
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this is in contradistinction to the approach of generic P(X, <fi) models (DBI being clearly 
an exception) and to the effective theory of inflation in its most comprehensive disguise 
[29] . It is the combined (signatures & stability)-oriented approach which is at the heart 
of our work here. 

Folded configuration. The tetrahedron has k\ 2 = 0, k% = k 2 and k 3 = fc 4 . The 
variables are x = k±/k\ and y = k\±/k\ with domains [0,1] and [0,2] respectively (also 
1 — x < y < 1 + x). In this configuration, T\ produces a shape which is very similar to 
the ones we can observe from contact interactions in P(X, <fi) models, whereas T 2 is a 
slightly different version of it in such that, unlike J-\, it is not constant along the x = 1 
axes. T 2 also shares some features with the shapes from scalar-exchange diagrams in 
P(X, <fi) models (also characterized by a non constant values for x — 1) however, unlike 
the latter, it is not null in (x, y) = (1, 0) or (1, 2). 

Specialized planar configuration. In this limit ki = k 3 = /c 14 and k± 2 can be expressed 
in terms of x = k 2 jk\ and y = k^/kx (both defined in the interval [0, 2]) 



We considered, without loss of generality, the positive sign solution for our plots. Notice 
that, again, T\ is very close to the contact interaction shapes for P(X, <ft) models, 
whereas T 2 is quite similar to the shape functions arising in these models from scalar- 
exchange diagrams. 

Planar limit double- squeezed configuration. In this configuration k^ = k^ = k\ 2 and the 
tetrahedron is going to be flattened. This allows to write k 2 in terms of x = k 3 /ki and 
y = ku/h 



where x G [0,1], y £ [0,2] and 1 — x < y < 1 + x apply. T\ has a shape that is very 
similar to the one from contact-interaction diagrams in P(X, <fi) models and it becomes 
null as x — > 0; the same behaviour in the x — > limit is observed for T 2 \ however, 
the latter slighlty differs in such that, unlike the former, it goes from being positive to 
being begative along the x — 1 axes, a feature which is new with respect to the existing 
literature. 

To summarize our findings about shape functions, T\ noticeably reproduces the features 
that can be also observed in the shapes from contact-interaction diagrams in P(X, <ft) 
models (and, thus, in DBI). 

As far as T 2 is concerned, even though it arises (like T\) from contact-interaction 
diagrams, it presents similar features to the scalar-exchange diagrams in P(X, 4>) models 
in the planar configuration and, but only to some extent, also in the folded configuration 




(53) 



1 + x 2 + y 2 + VS^J-I + 2x 2 + 2y 2 + 2x 2 y 2 - x 4 - y 4 , (54) 



14 



1.0 



Figure 2. Plot of T\ in the equilateral (upper left), folded (upper right), specialized 
planar (lower left) configurations and plot of Fi / {kik 2 k 3 k 4: ) in the planar-limit double 
squeezed configuration (lower right). 



and in the planar limit of the double- squeezed configuration (in the latter two, T 2 
partially differentiates itself from the P(X, <fi) results in some regions of the domain). 
Intriguingly, the T% shape in the equilateral configuration cannot be obtained at all in 
the P(X, <f>) context and therefore represents a signature for our type of models, quite 
a unique one in the realm of stable theories. 

Let us now compute the trispectrum amplitude. Following the convention in [17], 
we can quantify it by considering the value tnl m the limit (in (ITJ dubbed as RT, 
regular tetrahedron, limit) where k\ = k 2 = k 3 = k A = fc 12 = k u = k || 



where is the amplitude of the power spectrum of primordial scalar perturbations, 
which in our case is = (H 2 )/(eMp) . The three independent contributions to the 
total value of t^l read 




(55) 



r (i) 

1 NL 



189 H 2 Q 1 
512 eM| ' 



r (2) 
' NL 



4815 H 2 Q 2 
8192 eM 2 P ' 




1509 H 2 Q A 
16384 eMj, ' 



(56) 



§ Notice that t^l corresponds to the £nl quantity of [T7] . 
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For detection in the CMB, t^l needs to be several orders of magnitude larger than /jvl- 
If we consider our Lagrangian (e.g. as in Eq. (25)), we can see that 

f NL ( ~ g ~ (!) , r NL e ~ I ~ (f ) , (57) 

so the bispectrum is naturally easier to detect than the trispectrum. However, as the 
authors of [2T] observed, if we compare our setup to theories with dimension six operators 
at the level of the cubic Lagrangian, dimension eight for the quartic Lagrangian, e.g. 
P(X,<j>) inflation HUES], we have 

f NL c ~ 1 ~ (f ) , r NL e ~ 1 ~ (f ) . (58) 



Putting Eqs. (57 58) into words: for any given value of /jvl, a larger tjvl would arise in 
our model; precisely we have 

/ f 4 \ 2 / 5 

TNL~ ^ , (59) 



where we set ( 2 ~ ^ ~ 2 x 10~ 9 . The message is that this model is generally favored 
in comparison with models described by lower dimension operators (DBI being one of 
them), when it comes to the possibility of detecting tjvl. [jj 

As an order of magnitude estimate, for a bispectrum non-linearity parameter |/nl| 
of order 1, 10, 50, and 100 one expects a \t nl \ trispectrum amplitude of the order 50, 
2 x 10 3 , 3 x 10 4 , and 9 x 10 4 , respectively, g] 



|| For constraints on some specific forms of primordial trispectra, see, e.g. 1571 [551 15], 

% The values of /nl chosen for this simple estimate are within the 95% C.L. constraints set by Planck 

on the equilateral and flat bispectra, which are those generated by the class of models studied here [3T] . 
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Figure 3. Plot of T% in the equilateral (upper left), folded (upper right), specialized 
planar (lower left) configurations and plot oi T2/ {kik^kzk^j in the planar-limit double 
squeezed configuration (lower right). Notably, the equilateral configuration plot is 
unlike any other obtained within P(X, <f>) models. 

4. Conclusions 

We have studied higher-order curvature correlators of a very general Galileon model 
of inflation. Despite its generality, this model retains all the essential features of 
Galileon models: it is stable in a quantum field theory sense and it is predictive in 
that interactions coefficients do not receive large corrections from renormalization. 

The key difference with respect to a seemingly analogous [TTl |2"U] realization 
is that the model investigated here is better described as a "Galileon theory of 
fluctuations' [21], as opposed to a Galileon theory of the full action. The non- 
renormalizations properties (with additional assumptions [21]) are expected to hold 
also in this model of perturbations around an FRW background. 

The Lagrangian comprises more operators than that of [13 [20], allowing for generic 
(y 2 ip) n interactions, of which the original Galileons of [19] are a small subset. As 
a consequence, the equations of motion are no longer immediately second-order in 
time derivatives. A field redefinition restores second order e.o.m. order by order in 
perturbation theory. This translates into the theory not being valid above a scale A 
(there would be ghosts), which is a perfectly acceptable fact considering that we are 
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seeking an effective field theory valid in a specific regime, around a specific background. 
On the other hand, having allowed for more general interactions, one is rewarded with 
a richer spectrum of non-Gaussianites. 

The model displays interesting non-Gaussianities already at the level of the 
bispectrum [2T], but, considering that other more or less stable models do the same, 
it is important to extend the analysis to the trispectrum so as to remove as many 
degeneracies as possible in the predictions. Doing so does indeed pay off: in the 
equilateral configuration one finds a shape-function which cannot be generated in the 
entire class of P(X, <p) models, of which DBI-inflation is the most notable example. 

As for the impact of Planck data [5] on primordial non-Gaussianities, the simplest 
single-field models of slow-roll inflation are in agreement with observations, thus passing 
the most stringent tests of Gaussianity performed to date. However, for what pertains 
alternative signatures, a window for /tvl is still allowed, where it might be small 
but not zero and might peak in configurations which are not typical of the simplest 
inflationary models. Also, if one is after stable models of inflation, then the number 
of available options shrinks considerably and the one described here is certainly one of 
them. In light of the very recent results of the Planck mission [5J, we stress here that 
the bispectrum predictions of the models studied here comprise an equilateral as well as 
a flat contribution, which are consistent with the constraints provided by the data [5]. 

As mentioned, tnl detectability is also slightly enhanced in this model with respect 
to P(X,4>). This is especially interesting if one is willing to consider the possibility of 
a naturally small three-point function combined with a large four-point correlator [54J . 
Although this does not occur in our specific model (we restrain from fine-tuning the free 
M n coefficients), it is nevertheless a feature most likely to be realized in models highly 
constrained by symmetries, such as the one studied here. 

As for increasing the sensitivity on measurements for non-Gaussian observables, 
one expects the inclusion of polarization (and more data) in Planck data analysis to 
improve existent constraints. The same expectations hold true for LSS, see e.g., [59| 160] 
and possibly with CMBPol- [SI] and Core-like jH2] future missions, as well as 21cm 
cosmology [631 EH E5] ■ 
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Appendix A. Field redefinition 



The total action (|7j) counts a quite large number of terms, therefore finding the right 7T( 2 ) 
requires a long calculation. One convenient way to proceed is to work order by order 
in the number of temporal derivatives and cure each order at a time, starting from the 
highest. For the highest order, which is with ten time derivatives, we have 

SFV(1)] = / dV {-^Mlml^f) , (60) 

sF°V(o^(i)] = / d 4 xa 3 (-6a) (M + M 5 ) (nf 0) (n^) 2 + . (61) 

SrVfo)]=0, (62) 



'(o)J 

where 7t|q| = d 3 ir^o)/dt 3 and tt|q] = d 4 n^o)/dt 4 . We used 7T(i) from (17), which can also 
be rewritten as 

tt (1) = a7T( ) + 3^7^(0)^(0) - 7*(o) ( ] > (63) 



where a = (M 1 + M 5 ) j2M 2 P E and 7 = (2Mi + M 2 + M 6 ) /2M%H). After integrating by 
parts and summing up the above contributions, the fourth order action with ten time 
derivatives (iS>4 2 '[7r( )7r( 2 )] excluded) reads 

fi f = / (ftM .«).(5^*£). (64, 

These contributions are all unwanted and can be canceled by some of the terms 
[ 7r (o) 7r (2)] by choosing 



_ (3) , 2 / 5(Afi + M B ) 



2 % 



2M^7r (2) D -n (0) (^>Y -^w^ ) • (65) 



Mf>H 



This first contribution to 7T( 2 ) will also produce contributions to <S| and «S| from 

<S 2 2 V(o)] D J d A xa 3 2M 2 P Hv {2) (^Hir {0) - . (66) 

The next step will be to write down the complete expression for Sf and, similarly for 
what we just did for Sf°, find the contributions to 7T( 2 ) that will cancel its unsafe terms. 
We can see that 7i7 2 ) will receive new contributions at each order in t n and that, if one 
proceeds as for t 10 , the lower order will not affect the upper ones. The total expression 
for 7T(2) then counts several contributions that can be ordered as follows 

7T(2) = K(2,t™) + 7T(2,i 9 ) + 7T( 2ji 8) + 7T(2,t 7 ) + ^(2,t 6 ) + 7T( 2)t 5) + 7T( 2> t*) + ^(2,fi), (67) 

where 

2M 2 P HTT {2:tW) = n (0) (*§]) 2 C t u>, (68) 
2M 2 P Hn (2 ^ = 7r (0) (7rgj) 2 C t 9, (69) 

2M 2 p Htt m = 7ff 0) C ljt8 + Tr^TrgjC^ + (vrgj) 2 ^)^,^ + 7r (0 ) ^^ <V, (70) 
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9A/f 2rV _ -.2 . n , ■ _(3) dfo(O) r< , . (^Tf(O)) 2 ^ 

2MpHn(2,tf) = 7r(o) 7r (o)^i,f + ^(0)^(0) 2~^ 2 P + 2 Gs >* 7 

(X (X 
}2„ 



... ((9^ (0) )(^7r (0) ) (4) d 2 7T(o) r 

+7r (°) 12 64 >* 7 + ^(0)^(0)— ^— C 5,* 7 > ( 71 ) 

9A/f 2iV _ .2 r , -2 ^(Q) ^ , • « jjfofo) r , «2 ^(0) r 

ZMpUTT^t 6 ) = 7r(0)7T(0)^l,t 6 + ^(o) ~ ° 2 >* 6 + ^(0)^(0) — 6 3jt 6 + 7T(o) 2~~ G4 -* 6 

^(0)^(0) / gj^p) \ 2 / ^7T( ) \ / ^(Q) \ 2 

+vr (0 ) ~ 2 C 5>t e + vr {0) ^— j C 6 , t e + ^— j ( — j C 7 , t e 

+-(0) J J <V + 7T (0) ^ — J 6 9 , t e 

+7r (0) I Q 2 J I Q 2 J C 10.* 6 + >) I Q 2 J + ^(0) I Q 2 ) C U,*> 



(72) 



9/lf 2,v _ • ■■ dfo(o) r , -2 d 2 ^(Q) ^7r (0) gjjQ) 

2M p H 7r i2jt 5 } = 7T (0 )7r (0) — j—Ci^ + 7T( )— 2^ C 2,^ + 7T (0 ) — ^ ^ C 3,t 5 



cr w a" cr cr 

djS^difi^) <9i7r(o)^7f ( o)a 2 7T(o) ^^(o)^ 2 ^(o) ^ 
+7r (o) -4 6 4i t5 + — 6 5it 5 + 7T( ) — A 6 6jt5 

Ci Ct Lt 

+^(0)^(0) — -4 — 6 7;t B + 7T (0 ) -4 6 8jt5 + — 6 9jt 5 

Lb Lb Ct 

+7r (0) ^ <W' + 7T ( o) I I 6 11;t B + 7T (0 ) I q2 I 6i2, t 5, 



(73) 



ZMpHlTQt*) = 7T( ) 3— 6i )t 4 + 7T( ) ^ 6 2)t 4 + 7T( ) ^ 6 3jt 4 

a \ a § Cb 



n (0) —4 C 4,t 4 + ^(0) — 4 <-5,i 4 + ~ 4 6 6;t 4 

(V(o)) 2 ^7r ( o) 9i7r(o)^7r (0 )9i9j7r ( o) ^(0) (^(0)) 



^7,t A H g ^S,* 4 H 6 ^.t 4 



a 4 '' r a 6 a 6 

... df7i {0) dp 2 k 7r {0) di^didj^djdl-n^) 
+7T(o) ^ <~ao,t 4 H ^ °n,i 4 

^7T( ) (didj-kp)) d j n {td) d k d rn 'K( S)) d j d k d m TT iS)) 
H r W5,t 4 H p Oi 6 t 4 

| d^ {Q) djdjT: {0) did^ {Q) c ^^ ^ 

Cb^ 

2M p H7T(2,t 3 ) = - Q 6 lit3 + — 6 2 , t 3 + - 6 6 3i< 3 
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djit^didjitf^didl-Kio) dl^io) (pidj-K^ d j 7rf )d k di'Kr 0) d k d i d j 'Kt 0) 

+ ~ 6 C 4)t 3 + C 5>t S + C 6 ,3, 

(75) 

where the coefficients (T^n are linear combinations of H n Mi and H m Mij, where My = 
MiM^/MlH. 

We have thus shown how to perform the field redefinition up to fourth order. If we 
neglect all operators that are proportional to M (see discussion in Sec. 2.2.2), 7i7 2 ) reads 

- 2M 2 P Hn {2) = C t snf 0) + (C t r -3C t s) Htt^ + (C 1)t e + C#) tt 2 ' ^ W 



a 2 



+ (C 2 ^ + 3 (3(7(8 - C t i)) if 2 7r ( 2 0) 7f (0) + (C f 5 - 6C(8 + C t 7 - 3(71^6^/7 7T( )7f (0) 

(didjir^y (<9f7r (0) ) 2 
+ Ci,i 4 ^(o) "i 1 - (^.t 4 + Ci,t 6 + (7ts) 7f (0) — 

+ ^3,t 4 -" 7T(o) 2 r Oi )t 3.H7r(o) j h 2 ,t3-H 7T( ) j 

where 



C t8 = Mm + M n + M 3 , 
C t 7 = 3 (2Mn + Mi 3 + 4M 3 + M 7 y ■ 
C 1)t 6 = - (2M n + M 13 + 4M 3 + M iy 
C 2)t 6 = 3 (-2M 10 + 2M n + 3M 12 + 3M 13 + M 14 + M 2 + 18M 3 + 3M 6 + 9M 7 

C i5 = 2 (2M 10 - 2Mn - 3M i2 - 3M i3 - M 14 - M 2 - 18M 3 - 3M 6 - 9M- 
Ci )t 4 = -2M W - Mn + M14 + M 2 , 
C 2it 4 = M n + M12 + M X3 + 6M 3 + M 6 + 3M 7 , 



(77) 
(78) 

(79) 

(80) 

(81) 
(82) 



C 3it 4 = -9 (M 9 - M 10 - Mn 



3M 
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Ci it3 ee 3 (M 9 + 2Mi - Mn - M13 - M w + M 2 ) , 

C 2it 3 ee 2M 9 - 5Mi - 2Mn - 9Mi 2 - 2M i3 - 2M i4 + 2M 2 + 9M 3 + 9M 5 + 9M ( 



<7i,t 2 



-Mg - 2Mi + Mn + M13 + M14 - 

Mio + M12 - M 3 - M 5 - M 6 - M 7 . 



(83) 

Mi 3 - M M + M 2 + 3M 3 + 3M 5 + 3M 8 + 3M 7 ) , 

(84) 
(85) 
+ 9M 7 , 
(86) 
(87) 
(88) 



M 2 , 



which leads to Eq. (28). The coefficients Qi in Eq. (28) are 
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Qi = Mio + M12 + M 4 , 
Q 2 = M 8 - 2M 10 - M 12 + M 14 , 
g 3 = - 2M 8 - 8M 10 - IOM12 - 2M 14 - 12M 4 , 
Q 4 = - 4Mi - 6M 8 + 8M10 + 6M 12 - 2M ; 
Q 5 = 4Mi + 15M 8 + 28 



5Mi + 39Mi 2 + IIM14 + 54M 4 , 
g 6 = 6M1 + 9M 8 - 12M 10 - 9Mi 2 + 3M i4 , 



g 7 = - 12Mi - 36M 8 - 9M 9 - 39Mi + 9M n ~ 45M i2 + 9M i3 

- 27M 3 - 108M 4 - 27M 5 - 27M 6 - 27M 7 , 
Q 8 = M l + Mio - Mi 4 , 

g 9 = + 27A? 8 + — - -i 



15Mi 4 - 9M 2 



io 81M n 27M : 

* _~ 4 4 
243M 3 - 243M 5 243M 6 243M 7 

+ — + 81M 4 + — + — + — r^- 
4 4 4 4 



(89) 
(90) 
(91) 
(92) 
(93) 
(94) 

(95) 
(96) 



: 9 63M : 



81M 



4 



9M 14 81M 2 



(97) 



Appendix B. Computation of the Trispectrum diagrams. 

In the Schwinger-Keldysh (also dubbed as in-in) formalism, the general formula for the 
expectation value of a cosmological observable at a given time t is given by 



(©(*)> 



T(e i ti Hl{t ' )dt ' 



T e 



-i J*Hj(t")dt' 



(9£ 



where T and T are time-ordering and anti-time ordering operators and Hj is the 
interaction Hamiltonian (the subscript I indicates that the fields are in the interaction 
picture). Eq. (|98|) can equivalently be rewritten in terms of + and — fields 



<©(*)> = (T 



(99) 



In particular, we want to compute the expectation value at late times for trispectrum 
of the curvature fluctuations 



(100) 



(kkkk) = ( 27r ) 35(3) (E T c ^ fe, h) e(t) = 4444- 

Our interaction Hamiltonian is a function of the field operator rr (and its derivatives), 
which is related to the curvature fluctuation by ( = —Hit. We Fourier expand 

= I 7^**ftiE(*)> ( 101 ) 



(2vr)= 



where 



%(t)=ag7r fe (t) + a^(t) (102) 
and the creation and annihilation operators obey the standard commutation relations 



a h a k' 



(2tt) 3 ^ {%-%'). 



(103) 
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The Feynman rules for contracting the fields read 



T 
T 
T 
T 



7l + (x, t)lT + (x , t') 
7T + (x, t)7T~(x , t') 
7r~(x, t)ll + (x , t') 
7T~(X, t)-K~ (x , t') 



-> (n(x,t)7r(x,t'))9(t-t') + (n{x,t')n(x,t))6(t'-t), (104) 
-> (7r(x,t')7r(x,t)), (105) 
-> (ir(x,t)ir(x,lf)), (106) 
-> (ir(x,t)ir{x,t')}0(t'-t) + {Tr(x,t')Tr{x,t))6(t-t'), (107) 
and the external lines are to be treated as + fields. 

The connected trispectrum to tree-level requires the expansion of the exponential in 



Eq. (99) up to second order 



-i f*dt'(H+(t>)-Hy(t>)) 



+ 



J q dt {Hf{t')-Hj{t')) (10 
J^dt (Hf(t f ) - H Y (t')) J*dt" (H+(f) - Hj{t" 



2 Jo s ' ' ■> Jo 

and of the interaction Hamiltonian up to fourth order, i.e. 

Hj = J d 3 xa 3 (n {3) 



+ Hint) 



(109) 



(3) 

with two vertices, from T-L int . We will refer to these as contact-interaction and scalar- 
exchange diagrams. The latter turn out to be subdominant compared to the former (see 
discussion in Sec. 3), therefore we focus on the contact-interaction ones. Their formal 



Two different diagrams arise from (|108|), respectively with one vertex, from Hj^ t , and 



expression is as follows (we pick sample interaction term Ti^t = MH 4 tt 



4~4\ 



T™- (k u k 2} k 3} fc 4 ) ~ 2H 4 lm 



7r kl (t)7r k2 (t)7r k3 (t)7r k4 (t) J q ^VM#X(0^ 2 (*X 3 (0<(i 



Our eigenfunctions (in conformal time) are 

[1 + ikr) e-' ikT 
2y/ek 3 ~M P 



7Tfc( T ) 



111^ 



The time integral in Eq. (110) can be computed using the dimensionless variable y = k t r 



(where k t = k\ + k 2 + k 3 + k&) and then rescaling y with a parameter \i (this will be set 
equal to one at the end of the calculation) and noting that [21] 

i d 



TC. 



1 d 2 



■j-v-y 



• (113) 

and so on for higher powers of r. Note that the integrand functions are polynomials 
in kr multiplied by the phase e tktT . The polynomial part can be moved outside of the 



integral using (112) and (113). All we have left inside the integrals is 



dye 



1 



;ii4) 
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which we integrated using the standard prescription of performing an infinitesimal 
rotation of the contour along the imaginary axis at early times. The final results for the 
contact interaction diagrams from Eq. (28) read 



7> 01 = 



+ 



m*Q 1 

16fcifc 2 fc 3 fc 4 fc t 5 e 4 Mf, L" ' kf 
30 (fc 3 fc 4 + k 2 (h + fc 4 ) + fci (fc 2 + h + fc 4 )) 

T O 
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